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ABSTRACT 

We have studied the effects of electron-ion non-equipartition in the outer regions of relaxed clusters 
for a wide range of masses in the ACDM cosmology using one-dimensional hydrodynamic simulations. 
The effects of the non-adiabatic electron heating efficiency, /3, on the degree of non-equipartition are 
also studied. Using the gas fraction /gas = 0.17 (which is the upper limit for a cluster), we give a 
conservative lower limit of the non-equipartition effect on clusters. We have shown that for a cluster 
with a mass of Mvir ~ 1.2 x 10^^ Mq, electron and ion temperatures differ by less than a percent 
within the virial radius i?vir- The difference is ~ 20% for a non-adiabatic electron heating efficiency 
of /3 ~ 1/1800 to 0.5 at ^ 1.4i?vir- Beyond that radius, the non-equipartition effect depends rather 
strongly on /3, and such a strong dependence at the shock radius can be used to distinguish shock 
heating models or constrain the shock heating efficiency of electrons. With our simulations, we have 
also studied systematically the signatures of non-equipartition on X-ray and Sunyaev-Zel'dovich (SZ) 
observables. We have calculated the effect of non-equipartition on the projected temperature and 
X-ray surface brightness profiles using the MEKAL emission model. We found that the effect on 
the projected temperature profiles is larger than that on the deprojected (or physical) temperature 
profiles. The non-equipartition effect can introduce a ~ 10% bias in the projected temperature at 
i?vir for a wide range of /3. We also found that the effect of non-equipartition on the projected 
temperature profiles can be enhanced by increasing metallicity. In the low-energy band < 1 keV, the 
non-equipartition model surface brightness can be higher than that of the equipartition model in the 
cluster outer regions. Future X-ray observations extending to ~ i?vir or even close to the shock radius 
should be able to detect these non-equipartition signatures. For a given cluster, the difference between 
the SZ temperature decrements for the equipartition and the non-equipartition models, JATsze, is 
larger at a higher redshift. For the most massive clusters at z ~ 2, the differences can be SAT^ze ~ 4- 
5 /iK near the shock radius. We also found that for our model in the ACDM universe, the integrated 
SZ bias, Ynon-eq/Yeq, evolves slightly (at a percentage level) with redshift, which is in contrast to 
the self-similar model in the Einstein-de Sitter universe. This may introduce biases in cosmological 
studies using the /gas technique. We discussed briefiy whether the equipartition and non-equipartition 
models near the shock region can be distinguished by future radio observations with, for example, the 
Atacama Large Millimeter Array. 

Subject headings: cosmic microwave background — galaxies: clusters: general — hydrodynamics — 
intergalactic medium — shock waves — X-rays: galaxies: clusters 



1. INTRODUCTION 

Observational and theoretical studies have shown that 
the study of the intracluster medium (ICM) can be used 
as a test of plasma physics under extreme environments 
that cannot be achieved in terrestrial laboratories, as well 
as an important cosmological probe. If we assume the 
matter content of clusters is a fair sample of the universe, 
the baryon fraction of clusters can be used as an esti- 
mator of the average value for the universe, with proper 
correction for the baryons contained in the stellar compo- 
nent and for a small amount of baryonic matter expelled 
from clusters during th e formation process. Recent work 
from ( lAllen et al.ll2008[ ) has shown that the combined re- 
sults of the baryon fraction from X-ray observation with 
the cosmic microwave background (CMB) data can give 
powerful constraints on cosmological parameters, such as 
the equation of state of the dark energy. 

However, the study of cosmology using clusters of 
galaxies relies heavily on the understanding of cluster 
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physics. For precision cosmology, systematic uncertain- 
ties at even the percent level are significant. For example, 
in current studies, the baryon fraction within clusters is 
assumed to be independent of redshift and the mass of 
clusters. It would be important to see if these assump- 
tions are justified; if not, it is important to study the 
dependence of the baryon fraction on cluster properties 
and redshift. Even if the dependence on redshift is weak, 
the correction factor for the baryon content within clus- 
ters compared to the average value in the universe could 
affect the constraints on cosmological parameters. 

Studying cluster outskirts (> i?20(i_l) is very impor- 
tant because the boundary conditions of the cluster out- 
skirts constrain t he global properties of a cluster (e.g., 
iTozzi et all l2000f ). Also, the outer envelopes of clus- 
ters have been thought to be less subject to some ad- 

^ -Ra is the radius within which the mean total mass density of 
the cluster is A times the critical density. The virial radius -Rvir is 
defined as a radius within which the cluster is virialized. For the 
Einstein-de Sitter universe, -Rvir ~ Rns, while for the standard 
ACDM universe, Rvir ~ -^^95- 
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ditional physics including active galactic nucleus (AGN) 
feedback, and that the outer regions of clusters may pro- 
vide better cosmological probes. Currently, there are 
very few observations of the properties of the ICM in 
the outer parts of clusters. Thus, most of our under- 
standing of these regions is still based on numerical hy- 
drodynamic simulations which assume the hot plasma 
is a fluid. In these simulations, the clusters are formed 
from mergers and accretion of dark matter and baryonic 
gas in overdense regions. A variety of shocks with differ- 
ent geometries along the large-scale structure (LSS) fila- 
ments and transverse to them near and beyond the virial 
radius are unambiguous predictions of the cosmological 
hydrodynamic simulations. Unfortunately, the lack of 
observational information on the clusters outskirts pre- 
vents us from understanding the accretion shock region, 
and hence the input physics for the numerical simula- 
tions is called into question. For example, the thermo- 
dynamic state of the shocked g well as the shock 
position, depends on the pre-shock gas temperature; the 
shock will be wea ker if the infalling gas is pre-heated 
(|Tozzi et al.ll20 06h. Even worse, recently it was noted 
that the non-fluid properties may be important in re- 
gions near the virial radius, where the Coulomb colli- 
sional mean fre e path is co mparable to the cluster size 
of a few Mpc (jLoebl 120071) . The Coulomb collisional 
timescale can also be of the order of the age of the cluster. 
This suggests that a full kinetic gas theory is needed in- 
stead of the fluid approximation when studying the gas 
properties near the edge of the cluster. Direct conse- 
quen ces include non-equipartition between elec trons and 
ions (|Fox &: Loe blll99 7l: lEttori fc Fabianlll998D . element 
sedimentation (Chuzho v fc Loebl |2004[ ) , and suprather- 
mal evaporation of hot gas from the clusters (|Loebll2007l : 
see also Medvcdcv 2007). Some of these effects can lead 
to a bias in baryon fraction measurements, and hence 
cosmological studies. 

Recently, progress has been made in the study the 
baryon content of the outer regions of clusters through 
the X-ray observations together with the Sunyaev- 
Zel'dovich (SZ) effect on th e CMB by the hot e lectrons 
in the ICM out to - i?2oo (|Afshordi et al.ll2007| ). A 3ct 
result from the WMAP three-year data suggests that 
35% ± 8% of the thermal energy in ICM are missing, 
indicating that the baryons in clusters may be missing 
even accounting for those locked in stars. The result is 
also supported by independent measurements from other 
X-ray and SZ observations (Ettori 2003; LaR ogue et al.l 
[200a iVikhlinin et"all [200l lEvrard et ah .2003). Using 
the X-r ay observatio n s tog ether with numerical simu- 
lations, lEvra rd et al.l (|2008[ ) reported that as much as 
50% of th e thermal ene r gy can be missing in the ICM. 
Although iGiodini et al.l ()2009f ) reported that the total 
baryon fraction within a smaller radius of -R500 of mas- 
sive clusters are consistent to the cosmic value within la 
when all the X-ray hot gas, stellar mass in galaxies, gas 
depletion during cluster formation, and intracluster light 
from stars are taken into account, if the missing baryons 
measured in the outer region (< i?20o) is really signif- 
icant, this may indicate either a yet-unknown baryonic 
component, or some new astrophysical processes in the 
ICM which is driving out the gas from the clusters. While 
there is no evidence f or any undetected baryonic compo- 
nent, lAfshordi et al.l (|2007l) pointed out that the missing 



of hot baryons can either be explained by the thermal 
diffusion or the evaporation of baryons out of the virial 
radius of clusters. Another possibility is that electron 
temperature is low er than that of the cquipartition value 
dWon g et al.l[2008l) . 

Given the advancements in the X-ray and SZ obscrva- 
tions of the cluste r outer regions (lAfshordi ct al. 200^ 
Bautz et al.l 120091: iGeorge et all 12009^ .Reiprich et atl 
20091) ■ as well as the growing evidence of missing thermal 
energy in the ICM and the possible negative implications 
for cosmological tests, a more detailed study of the ki- 
netic processes in cluster envelopes is necessary. While 
magnetic fields may affect some of the kinetic effects of 
transport processes such as thermal conduction, the mag- 
netic effects on non-equipartition should not be impor- 
tant since the physics is local. Moreover, it is known 
that various astrophysical shocks in magn etized environ- 
ment lead to no n-equipartition (.Ghavamian et al.ll2Q07l : 
IHuh et al.l[200l . 

The collisionless accretion shock at the outer bound- 
ary of a cluster should primarily heat the ions since 
they carry most of the kinetic energy of the infalling 
gas. Assuming that cluster accretion shocks are similar 
to those in supernova remnants, the electron tempera- 
ture, Te, immediately behind the shock would be lower 
than the ion temperature, T^. The equilibration between 
electrons and ions would then proceed by Coulomb col- 
lisions. Near the virial radius, due to the low density, 
the Coulomb collisional timescale can be comparable to 
the age of the cluster, and the electrons a nd ions may 
not a chieve cquipartition in these regions (|Fox fc Loebl 
[T997I) . Since X-ray and SZ observations measure the 
properties of the electrons in the ICM, the net effect 
is to underestimate the total thermal energy content 
within clusters. This might account for some or all of 
the missing thermal energy in the ICM derived by the 
X-ray and SZ observations. As mentioned above, non- 
equipartition of ions and electrons is observed in var- 
ious astrophysical shocks. Most supernova remnants 
with high Mach numbers comparable to cluster accre- 
tion shocks have electron te mperatures which are lo wer 
than the ion temperatures (|Ghavamian m 
situ measurements from satellit es show the sam e fea- 
ture in the Earth's bow shock ()Hull et al.ll200H) . On 
the other hand. X-ray observations of the merger shock 
in the Bullet Cluster indicate that the equilibration time 
may be short er than that expected from C oulomb col- 
lisions alone (|Markevitch fc VikhlininI 120071 ). However, 
this merger shock has a Mach number of a few. From 
both s upernova remnant measu rements and a physical 
model, 'Ghava mian et ahl ()2007D have shown that elec- 
tron heating efficiency within a shock front (usually tens 
of the gyroradius) is inversely proportional to the Mach 
number squared. If the results can also be applied to 
the ICM, the low Mach number merger shocks would be 
immediately heated to T^/Ti ~ 1, while cosmological ac- 
cretion shocks with much higher Mach numbers would 
only be heated to T^/Ti ^ 1 by collisionless processes. 
After the electrons and ions pass through the thin shock 
front, they will likely be equilib rated by Coulomb colli- 
sions alone (|Bvkov et al.ll2008allH ). 

The non-equipartition in cluster of g alaxie s has 
been p reviously studied by iFox fc Loebl (|1997D and 
lEttori fc Fabian. (,1998. ) in semianalytic models. They 
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have shown that the temperature difference can be sig- 
nificant in the outer one-third of the shock radius of a 
cluster. One- or three-dimensional simulat ions for some 
indiv i dual clusters have also been studied (jChieze et alJ 
\199§; 'Takizawa 1999; Rudd & Nagai'2009'). While these 
simulations use different cluster or cosmological models, 
a general agreement is that the effect of non-equipartition 
is important if shock heating efficient of electrons is low 
1) and the equilibration afterward is due to Coulomb 
collisions alone. 

In this paper, we study systematically the effects of 
non-equipartition on X-ray and SZ observables in outer 
regions of relaxed galaxy clusters, which is particu- 
larly important for cosmological studies. We carry out 
one-dimensional hydrodynamic simulations with realis- 
tic Navarro-Frenk- White (NFW) model under the con- 
cordance ACDM cosmological background to provide a 
sample of clusters (groups) with different masses (10^^- 
lO^^M©) at different redshifts {z = 0-2). Even though 
we are studying the kinetic non-fluid properties in the 
cluster outer regions, the hydrodynamic treatment in 
modeling the cluster dynamical properties is reasonable 
and is justifled as follows. Even dynamically unimpor- 
tant magnetic fields should be able to reduce signifi- 
cantly the diffusion mean fre e path perpendicular to th e 
magnetic field (Sarazin 1984|Bor gam fc Kravtsovll200i ). 
The suppression of diffusion in a plasma depends on 
the topology of magnetic fields. For uniform magnetic 
fields, only diffusion perpendicular to the local magnetic 
field is suppressed, and along the field, particles move 
freely; their mean free path along a field line is still de- 
termined by Coulomb collisions. On large scales, diffu- 
sion is suppressed in bulk only if the magnetic fields are 
random and highly tangled on small scales. To include 
anisotropic diffusion in the calculation would be diffi- 
cult since the magnetic field structure is not known well 
enough. However, there is some evidence from large-scale 
magnetohydrodynamic simulations that magnetic fields 
in galaxy clusters are chaotic with correlation and rever- 
sal length scales o f ~ 50 and ~ 100 kpc, respectively 
(|Dolag et al.l 12002). Hence, we simply assume that dif- 
fusion is suppressed. We also assume that electrons and 
ions are equilibrated locally on a long Coulomb coUisional 
timescale, and assume that equilibration via plasma in- 
stabilities (Schckochihin et al. 2005, 2008) does not oc- 
cur except at the shocks (see Section[2]). Previous studies 
show that the dynamical properties of cluster outer re- 
gions in one-dimensional simulations successfully repro- 
duce those simulated in three-dim ensional calculations 
(jNavarro et al.l [T995t IRvu fc KangI 11997). The advan- 
tages of the one-dimensional simulations for our problem 
are presented in Section [5J We emphasize the signatures 
of non-equipartition on X-ray and SZ observations in our 
studies. We also study the effect of electron shock heat- 
ing efficiency on the degree of non-equipartition. Thus, 
observations of electron-ion equilibration may give con- 
straints to the electron heating efficiency, and hence the 
electron heating mechanism. The wider parameter space 
compared to previous work explored in this paper allows 
us to study the impact of non-equipartition effects on 
cosmological studies in a future paper. 

The paper is organized as follows. In Section [21 we 
describe the set up of our hydrodynamic models. The 



detailed implementations of the shock heating and the 
Coulomb equilibration process for our simulations are 
presented in Section |3l The ability of our simulations 
to reproduce analytic test models relevant to our stud- 
ies is discussed in Section ID We present the simulated 
dynamics of our realistic NFW cluster models in the 
standard ACDM cosmology in Section [5] These clus- 
ter models are used to study the non-equipartition ef- 
fects presented in the paper. We define the X-ray and 
SZ observables for our models to be studied, and also 
present the results for these observables in Section HI 
We discuss and conclude our work in Section [T] Un- 
less otherwise specified, we assume the Hubble constant 
Hq — 71.9 /i7i.9 km s~^ Mpc~^ with /171.9 = 1, the to- 
tal matter density parameter fluifl = 0.258, the dark 
energy density parameter fljy = 0.742, and the gas frac- 
tion /gas = flb/^M = 0.17, where flf, is the baryon den- 
sity parameter, for the realistic NFW model in the stan- 
dard ACDM cosmologjQ, and a hydrogen mass fraction 
X = 76% for the ICM throughout the paper. 

2. HYDRODYNAMIC MODEL 

LSS cosmological simulations predict that clusters do 
not evolve in isolation. During the linear phase of struc- 
ture growth, they can be influenced by tidal forces; while 
during the nonlinear growth phase, they can grow by ac- 
creting a signiflcant number of smaller clusters or merg- 
ing with clusters with similar sizes. However, here we are 
interested in studying the structure of accretion shocks 
in clusters, which can be observed most easily in re- 
laxed clusters which have not undergone a recent ma- 
jor merger. Cosmological studies using the gas fraction 
in clusters are restricted to clusters with the highest 
degree of dynamical relaxation to minimize systematic 
scatter in the dete rmination of cosmological parameters 
(| Allen et al.ll2008D : thus, relaxed clusters are of partic- 
ular interest. Moreover, based on a set of high reso- 
lution iV-body simulations, it has been found that the 
mass accretion history of a dark matter halo in gen- 
eral consists of two dis tinct phases: an early fast phase 
and a lat e slow phase (Wechsler et al.l 120021 : 1 Zhao et all 
I2003airbl : iLi et al.l 1200711 The fast accretion phase is 
dominated by major mergers, while the slow accretion 
phase is dominated by smooth accretion of background 
materials and many minor mergers. We are most in- 
terested in studying the non-equipartition of electrons 
and ions in the outer regions of clusters, where materi- 
als should be continuously accreting and the morphology 
is roughly spherical symmetric. Therefore, in our mod- 
els, we simply consider cluster growth by smooth accre- 
tion of materials from the background cosmology. In 
particular, we assume clusters are spherically symmetric 
and employ one-dimensional hydrodynamic simulations. 
It has been shown that one-dimensional calculations re- 
produce the density and temperature proflles of three- 
dimensional sim ulations of clusters in th e outer regions 
(jNavarro et al.r i995: R vu fc KangI Il997f ). where we are 
most interested. The assumption of spherical symmetry 
should be sufficient for us to gain insight into the astro- 
physical effects of non-equipartition on X-ray and SZ ob- 
servations, as well as the impact on cosmological studies 
(jWong et al.ll2008| ). Moreover, one-dimensional simula- 
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tions also allow us to better resolve shocked regions and 
to isolate the individual physical processes we are inter- 
ested in, so that non-fluid properties there can be studied 
in detail. This is difficult to achieve in three-dimensional 
simulations. Even though three-dimensional simulations 
have shown that accretion through filaments is a general 
feature in related clusters, these simulations also show 
that, other than the filament regions, material is accreted 
spherically and the morphology in the ou ter skirts of a 
clust ers are roughly spherical symmetric (|Molnar et al.l 
120091) . A very recent three-dimensional study has al- 
ready shown that the signature of non-equipartition for 
a relaxed cluster is roughly spherical symmetric as well 
(model CL104 in lRudd fc Nagail 120091) . 

2.1. Simulation Code 

We employed the PLUTO code (jMignone et al.l 120071 ) 

to solve the one-dimensional Newtonian hydrodynamic 
equations in spherical coordinates for our problem. This 
code provides a multi-physics, multi-algorithm modular 
environment which allows new physics to be included and 
new modules to be developed easily. The code does not 
include self-gravity or an A'^-body solver for determining 
the dark matter distribution, but these can be handled 
as force terms in the code. We have developed a scheme 
to include the dark matter contribution to the gravity by 
evolving the NFW profile sclf-consistently with the hy- 
drodynamic evolution of the fiuid calculated in the code 
(Section 12. 5p . This implementation has an advantage 
over using the A^-body solver, in that the dynamics of 
gas in the NFW dark matter potential can be investi- 
gated under controlled conditions. We have also imple- 
mented self-gravity of the gas in the code. Dark energy 
can also be included in the code easily as a force term 
(Section l2.6p . The code is built on modern Godunov-type 
shock-capturing schemes which are particular suitable for 
computation of highly supersonic astrophysical fiows in 
the presence of strong discontinuities. Shock-capturing 
is needed in our calculations since we are interested in 
calculating the non-equipartition signatures around the 
shock regions. 

2.2. Boundary Conditions and Computational Domain 

The boundary conditions depend on the geometry and 
the physics of the problem being solved. For a spheri- 
cally symmetric geometry of the ICM, the inner bound- 
ary condition is refiective. Because we are studying the 
smooth accretion of background materials in isolation, 
it is natural to use the Hubble-flow-like outer boundary 
conditions which are defined as 



dVg_ 




dlVg_ 


dr 


in 


dr 


Pa 


in — 


Pg, out 



and 

Pg.in Pg.out : (3) 

where r is the radius, Vg, pg, and Pg are the gas ve- 
locity, mass density, and pressure, respectively, and the 
subscripts "in" and "out" denote the inner and outer 
quantities at the outer boundary, respectively. Each sim- 
ulation is set up such that there is an over density in 



the central region compared to the critical density (Sec- 
tion 12. 3p . Near the over dense region, materials are ac- 
creted toward the center. To avoid boundary effects, the 
size of the computation domain should be large enough 
that materials should always be outgoing and follow the 
Hubble-fiow-like outer boundary condition. If the size 
of the boundary is too small, the dynamics of materials 
near the boundary would be influenced by the central 
cluster and would not follow the Hubble-flow-like out- 
flow; in the extreme case, materials should be infalling 
rather than outgoing. Bertschinger (1985) has shown 
that in an Einstein-de Sitter universe, materials follow 
the Hubble-flow-like outflow closely outside of a few times 
of the turnaround radius (Figure 1(a) therein). The 
turnaround radius is about 10 Mpc for a cluster with 
^ IO^^Mq at present, and hence the estimated radius 
at where materials follow the Hubble-fiow-like outfiow 
is about 30 Mpc from the cluster center. This condition 
should be sufficient for the ACDM universe, since the ma- 
terials in the ACDM universe should be less bounded. To 
be conservative, in our calculations, we take the bound- 
ary radius to be at least a few times this estimated radius. 
In particular, we take the boundary radius to be 100 Mpc 
in all models. The outer regions in our calculations at 
the final time reproduce the observed Hubble flow veloc- 
ity and the baryon density at redshift zero (Figures [T] and 
[5] in Sections |4] and m respectively), which validates our 
choice of the outer boundary conditions. 

The volume of each simulation was divided into 1000 
spherical annuli. The innermost zone width is set to be 
0.42 kpc. The widths of the zones increase with radius 
by a fixed ratio x, 

dr^+i = xdri , (4) 

where dri is the width of the ith zone. The grid is defined 
such that 

N 

Y^dn^L, (5) 

where A^ = 1000 is the total number of zones, and L = 
100 Mpc is the size of the computational domain. This 
is called a stretched grid in PLUTO. With this grid of 
radii, the zones near the accretion shock have widths of 
less than 1% of the shock radius in each simulation from 
z = — 2. Doubling the resolution gives essentially the 
same results. 

2.3. Initial Conditions for the Realistic NFW-Dark 
Energy Cluster Models 

Since at high enough redshift the dynamics of the back- 
ground universe is close to the Einstein-de Sitter model, 
the dynamics of a cluster should be close to a self-similar 
solution as well (jBertschingedll985l ). We chose a high ini- 
tial redshift and set up the initial conditions to be close 
to the self-similar solution. In particular, the initial con- 
dition is such that there are two regions separated by 
an accretion shock. The location of the cluster accretion 
shock in the baryonic material is very close to the first 
caustic in the dark matter, and the distribution and evo- 
lution of the baryonic gas and dark matter are almost 
identical outside th e accretion shock ()Bertschingeiill985l : 
iRvu fc Kandll99 7[). Within the shock radius, we assume 
the dark matter follows the NFW profile ((Navarro et ahl 
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[Till), 



Pdm{r) 



Pdn 



(r/r,)(l + r/r,)2 



(6) 



where pdm,s is a density scale and Vs is the scale radius. 
The scale radius, r^, is related to the concentration pa- 
rameter, c, and the virial radius, rvir, by 

(7) 



The virial radius, r^i 

Mvir 



., is defined by 
47r 

= — AvirPc(2)r: 



3 

vir ' 



(8) 



where Mvir is the virial mass, pdz) is the critical density 
of the universe at redshift z, and the critical overden- 
sity Avir is obtained from the solution to the top-hat 
spherical collapse mode l, and can be approximated by 
(|Brvan fc NormanlllQQSn 



IStt^ + 82x - 39x^ (for Qr = 0) 



(9) 



where x = [17m,o(1 + zf/E{z)^] - 1, E{z)^ = Om,o(1 + 
z)^ + nnil + zf rJA, and VLr = XjiH^Rf with R 
here equals the current radius of curvature of the uni- 
verse. We adopt a concentration parameter given by 
equation (|17p in Section[23J in which c is effectively equal 
to 4 in all the initial models considered. We also assume 
that the gas in the initial models is approximately in 
hydrostatic equilibrium with the dark matter potential. 
The choice of such initial conditions is convenient be- 
cause there exists an analytic solution for isothermal hot 
gas in hydros tatic equilibri um with the dark matter po- 
tential (Maki no et al.|[T998l ). but not the total potential. 
This is convenient for setting up our initial models. In 
the simulation runs, the self-gravity of gas is indeed in- 
cluded so that such initial models actually deviate from 
hydrostatic equilibrium slightly. We have checked that 
the numerical solutions we are interested in are rather 
insensitive to the set up of the initial conditions, as long 
as the initial models are approximately in hydrostatic 
equilibrium. Beyond the shock radius, both the dark 
matter and the gas follow the self-similar infalling solu- 
tions. The initial time (age of the universe) is chosen to 
be 29.4 Myr, which corresponds to a redshift of z = 70.6 
in the standard ACDM cosmology. The details of the 
setup are described below. 

We first define an initial cluster mass, Afsh,i, at the 
initial time chosen. The initial cluster mass is distributed 
within a n initial shock radi us, i?sh,i- In the self-similar 
solution (jBertschingei]ll985l ). i?sh i and Mgh i are related 
by 

-Rsh,i — Ascale''ta,i , (10) 

and 



47r 



(11) 



where rta.i is the initial turnaround radius, pc,i is the 
critical density at the initial time, and, Ascaic = 0.347 
and TOgcaie = 3.54 are the dimensionless scaled radius 
and mass in the self-similar solution at the shock radius, 
respectively. This implies that the initial average over- 
density within the accretion shock is 



Ash.i — Wscale/Agcale 



84.73. 



(12) 



Within the shock radius, we assume the dark matter dis- 
tribution is given by the NFW profile (equation [6]) . For 
a fixed value of concentration parameter, c = 4 assumed 
at high redshift (Section 12. 4p . the initial scale radius Vs.i 



can be solved by using equation (C9) in lHu fc Kravtso^ 
(2003.1 . which is 



i?: 



/.h. = ^/(l/c) 



(13) 



where /(x) = a;^[ln(l -I- x"^) - (1 + 2;)"^], x{J) in 
an accurate fit t ing f orm is given in Appendix C in 
IHu h Kravtsovl (|2003n . and Avir is given in equation ©. 
In general, c depends on mass, and can be calculated 
iteratively. The scale density pdm.s in equation (|6]) is then 
fixed by requiring that the total dark matter integrated 
to the shock radius is equal to (1 — fgas)Msh,i- 

Within the shock radius, we assume the gas is initially 
in approximately hydrostatic equilibrium (vg = 0) with 
the dark matter potential. Initially, the gas is assumed 
to be isothermal. The gas distribution can be solved 
analytically, and is given by (jMakino et al.lll998D 



p,(r)=p,,oe-^^^/Ml + 



276/(2r/rs 



(14) 



The central gas density pg^ and b are fixed by simulta- 
neously requiring that the gas fraction within the shock 
radius is equal to /gas and the gas density at the shock 
radius is given by the strong shock jump condition 



Pg2 = 4pgi 



(15) 



where the subscripts 1 and 2 denote the preshock and 
postshock quantities, respectively. Once b is fixed, the 
temperature, and hence the gas pressu re is indeed fixed 
by using equation (9) in Makin o et al.l ([1998). The gas 
pressure, Pg, is proportional to pg since the gas is as- 
sumed to be isothermal. However, if such pressure pro- 
file is used, the pressure at the shock radius does not 
match the strong shock jump condition in general. This 
is because the dynamical solution with accretion shock 
does not follow t he isothermal hydrostatic solution given 
m M akino et all ([T998). Since we are interested in the 
shock solution in our calculation, we renormalize the gas 
pressure such that the strong shock jump condition is 
satisfied at the shock radius. 



Pa2 = -Pgl Wgi 



(16) 



With this renormalization, the pressure is only lower 
than the isothermal hydrostatic equilibrium pressure by 
17%. The initial gas deviates slight from hydrostatic 
equilibrium. As long as the gas is in approximately equi- 
librium initially, the late time evolution of the cluster 
profiles in the outer regions we are interested in should 
not be affected by the slightly deviation in the initial 
condition. 

Gas in exact hydrostatic equilibrium with the total 
gravitational potential satisfying both the density and 
pressure shock jump conditions is not isothermal in gen- 
eral. We have checked that using initial conditions in 
exact hydrostatic equilibrium give essentially the same 
results. 
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Outside the shock radius, the gas density and velocity 
profiles are assumed to fo llow the self-similar infall so- 
lution (|BertschingeHll985l ). and the pressure is set to be 
effectively zero (i.e., the smallest value the code allows). 

2.4. Concentration Parameters 

We also need to set th e value of the concentration pa- 
rameter c. iBullock et al.l ()2001h proposed that c scales as 
(1 -I- z)~^. At low redshift, this is suppo rted by numeri- 
cal simulatio ns bv lWechsler et al.l (|2002| ) and lZhao et al.l 
(l2003bl l2009f ) , when the dark matter halos are in the slow 
accretion phase. However, at high enough redshift when 
the dark matter halos are in the fast accretion phase, 
the concentration parameter approaches a constant min- 
imum val ue independent of the mass and redsh ift of 
the halo (|Zhao et al.ll2003aL [2009t iGao et al.l[200l . At 
low redshifts, the variation of the concentration param- 
eter with m ass can simply be fitted using a power law 
()Dolag et a l. 2004; Neto ct al. 2007J Thus, we adapt 
the concentration parameter use d bvlBuo tc ct al.' (|2007[ ) 
but set a minimum value of 4 (jZhao et al. .2009.) inde- 
pendent of the redshift and the halo mass of the cluster 
in our simulations: 



c(Mvi,, z) 



Cm 

1 + 2 



M- 



14 



,4 



(17) 



where we take C14 ~ 9.0 and a — —0.172 as determined 
from the X-ray galaxy clusters with h alo masses between 
(0.06-20) X lO^^o at low redshifts (|Buote et al.l [20071) 

and Mi4 = 101-^/1^00^^0 • 

2.5. Self- consistent Evolution of NFW Potential 
Including Accreted Mass 

At each time step, we determine the gas shock ra- 
dius, and we assume dark matter is distributed accord- 
ing to the NFW profile within the shock radius; outside 
that radius, the dynamics of dark matter and gas are 
the same. We assume the baryon fraction is conserved 
as materials accretes within the shock, and hence, the 
dark matter mass within the shock radius is equal to 
(l//gas — 1) times the gas mass within that region. Al- 
though three-dimensional simulations suggest that there 
is a hydrodynamic outflow of gas and /gas inside a clus- 
ter is smaller than the cosmological value, our treatment 
here is at least self-consistent. Moreover, a smaller value 
of /gas inside a cluster in reality would enhance the ef- 
fect of non-equipartition, and hence our calculations give 
a conservative estimate of the non-equipartition effect. 
As mentioned in Section [^31 if the shock radius and the 
mass within it are known, the dark matter profile can 
be solved iteratively together with the concentration pa- 
rameter. The dark matter profile within the shock ra- 
dius calculated in each time step is used to calculate the 
gravity contributed by the dark matter within the shock 
radius. 

2.6. Dark Energy Implementation 

We will do the calculations in a cosmological back- 
ground including the accelerated expansion due to dark 
energy. The accretion of materials depends sensitively on 
dark energy since it governs the expansio n of the back- 
ground materials which are to be accreted (|Rvu fc Kand 



Il997f) . To implement dark energy in the one-dimensional 
Newtonian code, we utilize the similarity between the 
Friedmann solution and Newtonian cosmology, since the 
latter can be solved in the Newtonian code. The dy- 
namical solution of the scale factor, a, of a homogeneous 
universe is governed by 



47rG 



3P 



(18) 



where p and P are the total energy density and pressure 
of the universe, respectively, and G and c are the grav- 
itational constant and speed of light, respectively. The 
Newtonian analog is simply replacing a by r, i.e. 



df^ 



AttG 



3P 
72" 



(19) 



For a pure dark energy universe where the dark energy is 
described by the cosmological constant A, p ^ p\, P — 
Pa, and pA = — Pa/c^. Thus, equation becomes 



d r 
dt^ 



SttG 



-PAr . 



(20) 



We identify the right-hand side of equation (pOj) to be the 
repulsive force per unit mass contributed by the dark en- 
ergy. The dark energy density, pA, in the ACDM cosmol- 
ogy is constant throughout the history of the universe, 
and is given by 

PA = ^A,OPc,0, (21) 

where the quantities on the right-hand side are evaluated 
at the present time. With the force term in equation (|20p 
implemented into the code, the baryon density of the 
background cosmology is reproduced correctly in each of 
our simulations (see Figure [S] below) . 

3. ELECTRON HEATING WITHIN A CLUSTER 

3.1. General Picture of Shock Heating and the 
Equilibration between Electrons and Ions Afterward 

The heating of electrons passing through a shock 
within a cluster involves at least a two-step process. The 
first step is the shock heating of electrons within the very 
narrow shock, and the second step is the equilibration 
between electrons and ions afterward. In this section, 
we outline the general picture of electron heating within 
a cluster we assumed for our model. The detailed im- 
plementations of the shock heating (first step) and the 
Coulomb equilibration afterward (second step) are pre- 
sented in Sections 13.21 and 13. 3[ respectively. Section 13.41 
outlines the steps to calculate electron heating in our 
code. 

Immediately after the infalling material (electrons and 
ions, together with dark matter) has passed through a 
collisionless accretion shock, we assume it should primar- 
ily heat the ions rather than electrons since they carry 
most of the kinetic energy of the infalling gas. The elec- 
tron temperature immediately behind the shock would 
be lower than the ion temperature, and they would be in 
a state of non-equipartition. Indeed, non-equipartition 
of ions and electrons is known in various astrophysical 
shocks. Most supernova remnants with high Mach num- 
bers comparable to cluster accretion shocks have elec- 
tron temperatures which are lower than the ion tem- 
peratures (,Ghavamian et al..,2007i) ; in situ measurements 
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from s atellites show the same feature in the Earth's bow 
shock ()Hull et al.ll200lD . It is reasonable to assume that 
cluster accretion shocks are similar to those astrophysical 
shocks. Within the very narrow shock region, electrons 
and ions can exchange energy by coUisionless processes 
generated by plasma instabilities there. Unfortunately, 
the coUisionless heating mechanisms and the rate of heat- 
ing are still poorly known theoretically. We simply model 
the shock heating efhciency for the electron as a free pa- 
rameter (Section 13. 2p . 

After the electrons and ions have passed through 
the thin shock front , they will then be equilibrated by 
Coulomb collisions (jBvkov et al.l l2008al lbl). It is also 
possible that behind the shock, electrons and ions can 
exchange energy by coUisionless processes generated by 
plasma instabilities left behind. Again, the coUisionless 



heating mechanisms and the rate of equilibration are still 
poorly known theoretically. It has been argued that be- 
hind a strong shock, the spectra of supernova remnants 
are co nsistent with purely Coulomb equilibration pro- 
cesses (iR akowski et al.l l2003') . Given both the theoretical 
uncertainty concerning the coUisionless heating mecha- 
nisms and at least some observational evidence that col- 
lisionless processes are not important in the equilibration 
downstream, we simply consider the Coulomb collisions 
as the only equilibration process behind the shock in our 
calculations (Section 13. 3p . 

In calculating the equilibration processes of electrons 
and ions, we consider a fully ionized plasma of hydrogen 
and helium only. The timescale for charged particles of 
the same species to achieve a Maxwellian di stribution can 
be estimated by the collisional timescale (|Spitzeri[l962L 
p. 133) 



t^x = 28.6 Myr 



10 



r \ 3/2 , 

k) ( 



10- 



InA 



-1 ^1/2 



(22) 



where Tx is the temperature of the species x after it has 
achieved equilibrium, n,j. is the number density, is 
the particle mass number, Zx is the particle charge, and 
In A is the Coulomb logarithm which is similar for all 
species of interest. For the charged particles we are in- 
terested in, the electron-electron collisional timescale is 

1 /2 

the shortest, with tee tppAJ ~ ipp/43. For typical 
abundance nn ^ 10 riHe, tpp ^ iHc+2He+2. For pro- 
tons, the collisional timescale is tpp ^ 30 Myr, which 
is much shorter than the age of the cluster or the ac- 



I 

txy = 14.7 Myr ( ^^V^' f pL_ 

^ ^ yiO^Kj V10-5cm"3 



cretion timescale. For an accretion shock propagating 
outward with a velocity Vs ^ 1000 km s^^, protons at 
a distance d ~ Vstpp/A ~ 7 kpc away from the shock 
should be in a Maxwellian distribution. This is of the 
order of the zone size (^ 10 kpc) at the shock radius in 
our calculations. Hence, we can assume that protons, he- 
lium ions, and electrons are in equilibrium independently 
with temperatures Tp, Tfjg+2, and Tg, respectively. The 
energy exchange timescale for two species of particles, x 
and y in Maxwellian distributions with temperat ures 
and T y is given by the equilibration timescale (jSpitzeii 
[1961. p. 135) 

/lnA\"^ Ay ( TyA^y^'"^ 



Consider energy exchange between electrons and pro- 
tons, and assume that the electrons are initially much 
cooler than the protons. When Te < Tp{me/mp), the 
electrons are heated rapidly, with tep ~ ipp/3600. Thus, 
we expect that most of the time spent in electron heat- 
ing will occur in the regime in which Te/Tp 3> me/ nip. 
Hence, the energy exchange timescale between electrons 
and protons, tep is on the order of 0.5{mp/me)^^'^tpp ~ 
21tpp. Similarly, ^ tep for nn ~ 10 riHe- The 

equipartition timescale between electrons and ions (ei- 
ther protons or helium ions) will be of the order of 
0.6 Gyr, which is comparable to the accretion timescale 
of a cluster, [(cluster size)/ (accretion velocity)] ^ 1 Gyr. 
Thus, Te may not be equal to Tp or Tii^+2 in the outer 
regions of a cluster, which is the problem we are consider- 
ing. From equation (1^ . it can be seen that ^He+^p ^ tpp- 
It is also possible that Tp ^ 7hc+2 in the regions of in- 
terest. The solutions of the equilibration of the plasma 
formally can be obtained by solving the equilibration be- 



tween the three species. If we are interested only in the 
electron temperature evolution, we can consider the en- 
ergy exchange rate between the electrons and both of the 
two ion species 

l/te» = 1/tep + lAeHc+^ ■ (24) 

Since the contribution per ion to t~ ^ is the same for pro- 
tons and helium ions (equation [^3|) , the rate of change 
of electron temperature with the plasma, dTe/dt, de- 
pend s only on the mean ion temperature (|Fox fc Loebl 
|1997[ ). Moreover, once the ions reach equilibrium with 
each other, the rate of change of their temperatures due 
to collisions with electrons is the same, so they will re- 
main in equilibrium. Thus, for simplicity, we assume a 
single ion temperature, Ti, and we assume the electrons 
with temperature Te are equilibrated with this single ion 
temperature plasma. 

In hydrodynamic simulations, the three hydrodynamic 
variables, the gas mass density, pg, the total gas pres- 
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sure, Pg, and the gas velocity, Vg, determine the hydro- 
dynamic state of the fluid completely, independent of the 
kinetic state of the plasma. Since we are interested in 
the difference between electron temperature, Te, and io_n 
temperature, T^, we assign a new variable, r = Tg/T, 
in each grid of our hydrodynamic calculations, where T 
is the average thermodynamic temperature of the fluid 
given by T = {^mp/kB){Pg/ Pg), where /x = 0.59 is the 
mean molecular weight, rrip is the proton mass, and fcs is 
the Boltzmann constant. The electron temperature, Te, 
the ion temperature, T^, and the average thermodynamic 
temperature, T, are simply related by 



UeTe + UiTi 



(25) 



where rig and rii are the electron and ion number densi- 
ties, respectively. 

3.2. Electron Heating within the Thin Shock Front 

Within the very narrow shock front, electrons can be 
heated by both adiabatic compression and non-adiabatic 
heating. Non-adiabatic heating includes the conver- 
sion of the bulk kinetic energy into thermal energy and 
other collisionless heating processes generated by plasma 
instabilities. We can define the total non-adiabatic 
electron heating efficiency in the shock, f3, to be the 
change in electron temperature due to non-adiabatic 
heating, ATe^ non-ad, relative to the change in average 
thermodynamic temperature due to non-adiabatic heat- 
ing, ATnon-ad, 



/3 



V AT 



(26) 



non-ad 



The change in the bulk kinetic energy per particle is 



m(w2^ 



)/2, where m is the particle mass, and hence 



the increase in the electron temperature is much smaller 
than that of the ions in a shock. For the electron heating 
processes generated by plasma instabilities, the details of 
the mechanism and the efficiency are still unclear. There- 
fore, in our study, we assume total non-adiabatic electron 
heating efficiency, /3, to be a constant for any given clus- 
ter model. Adiabatic compression increases the electron 
temperature by a factor of {pe2/ Pei)'^~^ , where 7 is the 
adiabatic index, p^i and pe2 are the preshock and post- 
shock electron temperatures, respectively. For a strong 
shock with 7 = 5/3, adiabatic compression only increases 
the electron temperature by a factor of 2^/'^ sa 2.5. Thus, 
if the preshock electron temperature, Tei, is negligible 
compared to the postshock average temperature, T2, adi- 
abatic heating would not be important. This is indeed 
the case in our models. Nevertheless, we have included 
adiabatic heating in our calculations. Including both adi- 
abatic and non-adiabatic heating, the postshock electron 
temperature immediately after the shock is given by (Ap- 
pendix [X]) 



7-1 

Te2=[^] Tel + /3 max 



0,T2- ( ^ 



7-1 



(27) 

where a minimum of zero in the second term is set to 
ensure that numerical fluctuations do not introduce a 
false decrease in entropy in the non-adiabatic heating. 



Note that we assume that the heating of electrons is 
equivalent to the increase of temperature. This is only 
true if the electron distribution is Maxwellian. The time 
for the electrons to achieve a Maxwellian distribution af- 
ter they have been heated at the shock is tee ^ ^ 
0.7 Myr fSection 13. As shown in Section [3.11 this is 
much shorter than the accretion timescale, and the elec- 
trons within the zone size at the shock radius of our simu- 
lations should achieve the Maxwellian distribution. Note 
also that even if electrons are in non-Maxwellian distri- 
bution and if such electrons are heated much faster than 
the Coulomb coUisional heating rate, the faster heating 
within the shock has already been parameterized by the 
electron heating efficiency /3. 

Recently, observations of supernova remnants have 
shown that the electron heating efficiency at the shock 
is inversely proportiona l to the Mach number squared 
(jGhavamian et al.ll2007| ). If the result can be applied to 
cluster accretion shocks, the high accretion shock Mach 
number {M. > 100) would imply the electron heating effi- 
ciency to be <C 1 in cluster accretion shocks. In our work, 
we consider two cases for the shock heating of electrons: 
(3 = 1/1800 (the mass ratio of electrons to protons) as 
a model for a very low non-adiabatic heating efficiency 
which is supported by supernova remnant observations, 
and /3 = 0.5 as a model for an intermediate heating effi- 
ciency within the shock. 

3.3. Coulomb Equilibration After Shock Heating 

After the shock-heated material has passed through 
the thin shock front, we assume Coulomb equilibration as 
the only heating process for electrons behind the shock. 
The evolution of ele ctron temper ature due to Coulomb 
collision is given by (ISDitzedll962i p. 135) 



dr 
'dt 



2 In A 
503 



A 



T3/2 ' 



r-3/2n - 



(l-r)s- 



(28) 



where n is the total particle number density, InA ss 
37.8-Hln(Te/10SK)-ln(ne/10-3cm-3)i/2 is the Coulomb 
logarithm, the angle bracket term is the mean value over 
the ratio of the square of ion charge Z and the atomic 
number A. For our model of a pure hydrogen and he- 
lium gas, the angle bracket term equals 1 as long as 
Te/Ti > me/mp, which is true for our models with a 
minimum /3 = 1/1800. We solve equation (1^51) for each 
time step in our hydrodynamic calculations (SectionO 



3.4. Steps to Calculate the T^/T Evolution 

Cosmological accretion shocks are identified as shocks 
in which the Mach number, M, is greater than 10. In our 
simulations, we use the pressure shock jump condition. 



Pal 
P91 



7+1 7+1 



(29) 



to identify shocks. In each time step of our hydrody- 
namic simulations, if the condition of > 10 in equa- 
tion (|29l) is satisfied, we identify there is a strong shock 
at r2. We found that this condition can identify cosmo- 
logical shocks correctly in our simulations. _In regions 
without a strong shock, we calculate the Te/T evolution 
using equation (j28p throughout the entire time step. If 
an accretion shock is identified, we first use equation (|28p 
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to evolve T^/T for one half of the time step. Then, we 
apply the shock heating using equation (j27p . Finally, we 
use equation ()28|) again to evolve Te/T for the reminding 
half of the time step. 

Cosmological simulations show tha t there are variou s 
kinds of shocks formed in the LLSs ()Kang et al.l 120071 ) . 
Based on their location, these shocks can be classified 
as external and internal shocks. External shocks are 
formed around the outermost surfaces of the LLSs that 
unshocked intergalactic gas is falling onto (sheets, fila- 
ments, and halos), and the gas is shock heated for the 
first time. These external shocks in general have Mach 
numbers ^ 10. Internal shocks have low Mach numbers 
< 10. They are formed within those nonlinear structures, 
e.g., galaxy clusters and filaments, by the infall of pre- 
viously shocked gas during subclump mergers, as well as 
by chaotic flow motions. These cosmological simulations 
show that there are a large number of low Mach num- 
ber internal shocks contributing significantly to th e ther- 
mal energy budget of a cluster (jKang et al.ll2007[ ). Since 
those small Mach number internal shocks are mainly 
found well within inner regions of a cluster where the 
densities are high enough, electrons and ions are in gen- 
eral in equipartition. Here, we are only interested in the 
non-equipartition effects of electrons and ions, and these 
effects are only significant in very low density regions 
(e.g., cosmological accretion shock regions) rather than 
the contribution of the thermal energy from the shocks to 
the cluster. For our one-dimensional simulations which 
are used to study relaxed clusters, due to the symmetry 
of the problem, there is in fact only one cosmological ac- 
cretion shock with Mach number ^ 10 for each cluster, 
and there are no internal shocks formed. For compari- 
son, three-dimensional simulations have shown that for 
relaxed clusters, the number of small Mach number in- 
ternal shocks is much sm aller than that of th e high Mach 
number external shocks (jMolnar et al.|[2009l ). This is be- 
cause internal shocks are formed mostly by mergers or 
chaotic flow motions which are less likely to be found 
in relaxed clusters by definition. Thus, we do not have 
to consider electron heating by the low Mach number 
internal shocks at least for relaxed cluster that we are 
interested in. Future three-dimensional simulations will 
be needed to test our assumption. 

4. TEST MODELS: SELF-SIMILAR MODELS 

The essential physics governing the equilibration prob- 
lem we are considering are the hydrodynamics of the gas, 
the gravity by both dark matter and gas, the equilibra- 
tion physics between electrons and ions, as well as the 
dark energy which modifies the background cosmology, 
and hence the rate of accretion onto clusters and the du- 
ration of the accretion throughout cosmic history. The 
contribution of dark energy can be tested by whether the 
baryon density of the background cosmology can be re- 
produced correctly in each of our simulations, and this is 
indeed the case for our realistic simulations which will be 
presented in Section [5] (see Figure [5] below) . We present 
three tests here to show that our hydro code can handle 
the necessary physics correctly for the problem we are 
considering. The hydrodynamical response of the gas in 
an external gravitational field is tested by reproducing 
the analytic self-similar solutions of a collisionless dark 
matter dominated accretion model ([Bertschinger, ,1985i ) . 



The correct handling of self-gravity of the gas is tested 
by reproducing the analytic self-similar s olutions of a 
coUis ional gas dominated accretion model ([Bertschingeii 
|1985() . Finally, the correct handling of the equilibration 
physics is tested by comparing to the analytic solutions 
of a self-similar n on-equipartition model calculated by 
iFox fc Loebl (IT997I) . 

4.1. Collisionless Dark Matter Dominated Accretion 

Model 

In this section, we compare our numerical simulation 
with the analytic self-similar solution of the collisionless 
dark matter dominated accretion model in the Einstein- 
de Sitter universe {^m = 1,^^A = 0). In this model, we 
assume gravity is dominated by the dark matter poten- 
tial, and hence the self-gravity of the gas is switched 
off in the code. This is done to test the ability of 
our code to handle an external gravitational potential, 
which is important for our realistic simulations with the 
NEW dark matter potential included. We set fif, to be 
0.05 in this model. We assume the current Hubble con- 
stant to be i?o = 71.9 km s^^ Mpc^^ in this simulation 
run, and hence at redshift z — Q, the critical density is 
p^[z = 0) = 9.71 X 10~^° g cm~^. The initial gas pro- 
fil es are set according to the self-similar solution given 
in iBertschingeii ()1985[ ). The gravitational potential due 
to the dark matter is calculated according to the self- 
similar solution at each time step. The gas evolution 
within the dark matter gravitational potential is then 
calculated by the PLUTO code. Figure [1] shows the dy- 
namical variables of our numerical simulations for the 
self-similar model at four different redshifts. The model 
shown in the three panels is such that the dark matter 
mass accreted within i?i78 is Afiyg — 1O^^M0 at z = 0. 
The figure shows that a strong shock is propagating from 
i?sh ~ 0.7 Mpc at z = 2 to i?sh ~ 3 Mpc at z = 0. The 
shock velocity decreases from Vsh ~ 1500 km s^^ at z = 2 
down to Vsh ~ 1300 km s^^ at z = 0. At z = 0, the 
gas density just within the shock is about 0.8pc(2 — 0), 
and drops to about 0.2pc{z = 0) just beyond the shock, 
while a very sharp jump in pressure can be seen at the 
shock radius. At about 100 Mpc, which is at the edge 
of our simulation domain, the density drops according 
to the background cosmology as ^hPc (inset in the mid- 
dle panel of Figure [1]) . The velocity profile at very large 
radii also follows the Hubble fiow (inset in the top panel 
of Figure [1]) . 

To compare with the analytic self-similar solution di- 
rectly, we scale the physical radius and the dynamical 
variables to obtain the scaled dimensionl ess valuables 
according to the scaling relations given in IBertschingeii 
(1985,), 



A = 

m-- 



-r/rts, , 
_ Vgjr, t) 

_ Pg{r,t) 

Pa{r,t) 
Pc^brlJt^ 



(30) 



where A, v, 0, and ■0 are the dimensionless radius, ve- 
locity, density, and pressure of the gas, respectively, t 
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is the cosmic time which is equal to Xj yJ^ixGpc for the 
Einstein~de Sitter universe. The scaled dimensionless 
valuables in our simulation are shown as solid lines in 
Figure [21 The analytic solutions are plotted as open 
circles for comparison. After scaling, the dimensionless 
valuables in our simulation at the four different redshifts, 
as well as the analytic solutions, all lie almost along the 
same locus which can hardly be distinguished on the fig- 
ures. This shows that our numerical simulations are in 
excellent agreement with the analytic solutions within all 
region of our interest (~ 0.1-10 Mpc). We conclude here 
that our code can calculate the dynamics of the g as, as 
well as handling the shock in an external gravity under 
the cosmological expansion correctly. 

4.2. Collisional Gas Dominated Accretion Model 

Since we have included self-gravity of the gas in our 
realistic simulations, it is necessary to check whether 
our code can handle self-gravity correctly. We have per- 
formed a set of simulations to compare with the analytic 
self-similar solution of the collisional gas accretion model 
in the Einstein-dc Sitter universe {flM = = 1,^^A = 
0). In this model, the universe is considered as purely 
collisional fluid, and hence we do not include the dark 
matter contribution in these simulations. Self-gravity of 
the gas is included in the simulations. We scaled the 
physical radius and the dynamical variables according to 
equation ((30)) with fit, = 1. The scaled dimensionless 
valuables in our simulations are shown as the solid lines 
of Figure [31 The analytic solutions are plotted as open 
circles for comparison. Again, the figure shows that our 
numerical simulations are in excellent agreement with the 
analytic solutions. We conclude here that our code can 
calculate the dynamics of the g well as handle the 
shock with self-gravity under the cosmological expansion 
correctly. 

4.3. Self-similar Non-equipartition Model 

With the self-similar dynamical cluster model calcu- 
lated in Section |4 . 1 1 from our simulations, we can calcu- 
late the degree of non-equilibration of the gas following 
the method given in Section [31 In Figure [H the ratios 
of the electron and average thermodynamic temperature, 
T, for two clusters with different masses are shown as cir- 
cles. The correspo nding lines are the analytic solutions 
(|Fox fc Loebl 119971 ). Our simulations are in very good 
agreement with the analytic solutions. The slightly de- 
viations are mainly due to the finite resolution of the 
shock region in our numerical simulations. We conclude 
here that our code can calculate the equilibration physics 
correctly. 

5. DYNAMICS FOR REALISTIC NFW-DARK ENERGY 
MODELS 

Table [H lists the masses and radii for different over- 
densities for some representative NFW cluster models in 
the standard ACDM cosmology at z = 0. The different 
definitions of masses and radii are used interchangeably 
throughout the paper. In Figure [SI we show the evo- 
lution of the dynamical variables {vg , pg , and Pg) as a 
function of radius for a cluster with an accreted mass of 
Mvir = 1.19 X IQI^Mq at z = 0, where Mvir = M95 in the 
ACDM cosmology we assumed. The total mass within 
the shock radius is M^h = 1-53 x lO^^M© at z = 0. 



The figure shows that a strong shock propagates from 
Rsh ~ 1.2 Mpc at z = 2 to R^h ~ 4.2 Mpc at z = 0. The 
shock velocity decreases from Vgh ~ 1400 km s^^ at z = 2 
down to Vsh ~ 1000 km s^^ at z = 0. Similarly to the 
self-similar dark matter dominated model, at z = 0, the 
gas density just within the shock is about 0.8pc{z — 0), 
and drops to about 0.2pc{z — 0) just beyond the shock. 
At about 100 Mpc which is at the edge of our simulation 
domain, the density drops according to the background 
cosmology as fltPc (inset in the middle panel of FigurelS)). 
The pressure profile flattens as the NFW cluster evolves, 
and the central pressure drops. This is in contrast to 
the self-similar dark matter dominated model where the 
pressure never drops even in the very central region 
< 0.1 Mpc (not shown on the graph). The velocity proflle 
at large radius also follows the Hubble flow correctly (the 
inset of Figure [5]). This shows that with our implemen- 
tation of the dark energy, our simulations can reproduce 
the background cosmology correctly. For comparison, 
the shock radius and the shock velocity for a self-similar 
dark matter dominated cluster with the same mass (Mgh) 
in the Einstein-dc Sitter universe (J^m = 1,^^A = 0) 
are i?ss ^ 3.11 h^l(^^{MsY,/1.53 x IO^^Mq)^/^ Mpc and 
ufs ^ X 103 h^{%{Msb/l.53 X 1015Mq)1/3 km s^^, 
respectively. For the realistic NFW-DE model, the shock 
radius is larger than that of the self-similar solution by 
a factor of 1.38, while the shock velocity is lower by a 
factor of 1.35. The result is consistent with the one- 
dimen sional A^-body simulation given by iRvu fc Kand 
()1997|) . Their results show that the shock radius and the 
shock velocity for a cluster in the ACDM cosmology with 
r^A ~ 0.74 is about 1.4 larger and 1.35 lower than those 
in the Einstein-de Sitter universe, respectively. 

6. OBSERVABLES FOR REALISTIC NFW-DARK ENERGY 
MODELS 

6.1. Definition of X-ray Observables 

X-ray spectra depend mainly on electron temperature 
rather than ion temperature, and the projected tem- 
perature proflle of a cluster can be directly measured 
from X-ray observations. Different weighting schemes 
have been used to calculate the project ed temperature 
from hydrodynamic simulations, and iMazzotta et al.l 
()2004[) have shown that most of the commonly used 
weighting schemes (the mass- weighted and the emission- 
weighted) give signiflcantly different results from the X- 
ray observed spectroscopic temperature, Tspoc, and the 
spectroscopic-like temperature calculated from an ana- 
lytic weighting scheme they developed is able to approx- 
imate Tgpoc to better than a few percent for tempera- 
ture above ^ 3 kcV. However, for the non-equipartition 
model we considered, the electron temperature in the 
outer regions can be as low as 0.1 keV. X-ray line emis- 
sion from heavy elements contribute significantly to the 
X-ray spectrum for temperature be low ^ 3 keV . Such 
line emiss i on is n ot considered by Ma zzotta et al.l ([2004) . 
iVikhlininI (|2006D has generalized the spectroscopic-like 
temperature to lower temperatures down to ^ 0.5 keV 
and to arbitrary metallicity. This generalized scheme 
takes into account both the continuum an d line emis- 
sion a ssuming the MEKAL emission model ( Mewe et al.l 
[19851 [Kaaitra fc Mewel[l99l iLiedahl et al.l 119951) . The 
detector response has also been taken into account. The 
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weighting scheme is no longer analy tic and ha s to be tab- 
ulated. Following the algorithm of I VikhlininI (|2OO60 . we 
calculate the projected spectroscopic-like temperature, 
Tsi, from our numerical simulations. The integration is 
carried out within the shock radius. The latest Chandra 
ACIS-S aim point response files released for the ACIS Cy- 
cle 1 1 proposal planning are used to generate the weight- 
ing tableo Using different Chandra C CD response file s 
does not affect the results significantly (jVikhlininl I2Q06D . 

We also calculate the surface brightness profile from 
our numerical simulations. Here, we define x as the pro- 
jected radial distance from a cluster center. The surface 
brightness profile in a given energy band E is given by 



Ssix) = / AsiTe, Z)nenpdl , 



(31) 



where AE(Te,Z) is the cooling function which depends 
only on the electron temperature and heavy element 
abundances, Z, np is the proton number density, and 
/ is the distance along the line of sight. The integration 
is carried out within the shock radius. The cooling func- 
tion Ae(Tp,,Z) is calculated using the MEKAL model 
(iMewe et al.lll985t[Kaastra fc Mewill993l : iLiedahl et al.l 
[l995f K 

For our models with non-equipartition considered, T^. 



is used to calculate the T= 



sl,non-cq 



and Se. 



non-cq ■ 



For com- 



parison, we also consider models with electrons and ions 
are fully in equipartition by taking — T , which is usu- 
ally assumed in the literature. Quantities calculated with 
full equipartition assumed are denoted with the subscript 
"eq". 

Note that in the outer region where electrons and ions 
may not be in equipartition, non-equilibrium ionization 
may also be important. This may increase line emissions 
in the soft bands {E <\ keV). In our calculations, non- 
equilibrium ionization is not considered. 

6.2. Definition of SZ Observables 

The SZ effect by a cluster at x can be characterized as 
a temperature increment, ATszE(a^), with respect to the 
CMB spectrum 

ATsze(2:) = fi9)y{x)TcMB , (32) 

where 9 = /li^/fcsTcMB is the dimensionless frequency, 
TcMB is the CMB temperature, y is the Comptonization 
parameter, and f(6) is given by 



eO - 1 



(33) 



For simplicity, we neglect the relativistic term, ^sze, 
which is not important for frequency lower than about 
250 GHz. The Comptonization parameter is given by 



y 



nieC^ 



nJTedl (X / Pedl 



(34) 



where ctt is the Thomson scattering cross section and 
Pe = neksTe is the electron pressure. 

Similarly to the calculation for the spectroscopic-like 
temperature, for our models in non-equipartition and in 
fully equipartition, we calculate the Ar„on_(,q/cq(a;) by 

^ http://cxc.harvard.edu/caldb/prop_plan/imaging/index.html 



integrating equation p4p along the line of sight within 
the shock radius. 

Another useful SZ observable is the integrated Comp- 
tonization parameter, Y , which is defined as the integra- 
tion of the Comptonization parameter in equation p4p 
on the sky 



Y 



4 



ydQ 



ydA, 



(35) 



where dA is the angular diameter distance to the cluster, 
n is the solid angle of the cluster on the sky, and A is 
the projected surface area. We integrate the projected 
surface area of the cluster up to the shock radius. Such 
a quantity is useful for spatially unresolved clusters with 
SZ observations where the solid angle covers the whole 
cluster. 

Since Y = J ydA cx / PedV oc / n^TedV , where V is 
the volume of the cluster, the integrated Comptoniza- 
tion parameter Y is basically measuring the thermal en- 
ergy of the electrons in the cluster. If electrons and ions 
are not in equipartition and if the electron temperature 
is lower than that of the ions globally, the value of Y 
measured would be lower than the equipartition value. 
To characterize the degree of non-equipartition of the 
whole cluster, we define the bias as the ratio ynon-cq/^cq, 
where Ynon-cq is the integrated Comptonization param- 
eter for the non-equilibration model and l^oq is that for 
the equipartition model (Te = T). 

It has been shown that the integrated Comptoniza- 
tion p arameter displays a t ight correlation with cluster 
mass (|Reid fc Spergell [20061) ■ Such a tight correction is 
needed for precision cosmology, and hence a correct un- 
derstanding of the integrated Comptonization parameter 
is important. A detailed discussion of t he use of SZ sur- 
veys t o study cosmology can be found in lCarlstrom et all 
([200l . 

6.3. Results for the Temperature Profiles and the X-ray 
Observables 

The evolution of the average thermodynamic tempera- 
ture as a function of radius for a cluster with an accreted 
mass of Afsh = 1-53 x IO^^Mq at z = in the ACDM 
cosmology is shown in Figure [HI In general, for each red- 
shift, the temperature profile rises from the very central 
region to a peak, and then drops toward the outer region. 
The drop of temperature in the outer region is due to the 
drop of the shock velocity during the accretion history. 
The central drops in temperature, as well as the decrease 
of the value of the peak temperature are due to the adia- 
batic expansion of the cluster during the evolution. This 
is in contrast to the average thermodynamic tempera- 
ture of the self-similar solution in the Einstein-de Sitter 
universe where the temperature profile is always rising 
toward the central region at all redshift (Figure [T]) , and 
the temperature within the same gas mass, Mgas, never 
drops during the accretion history. The expansion of the 
hot gas in the central region of the cluster in our simu- 
lation is probably caused by the evolution of the NEW 
profile, since there is a central pressure drop for the gas 
in the NEW potential but not for the case in the self- 
similar model (Section [S]). In a real cluster, the temper- 
ature profile of the central region (< 100 kpc) is likely 
to be complicated by physical processes such as cooling, 
AGN heating (|Fabian et al.. .2000; .Blanton et al.. .20011 



12 



Wong & Sarazin 



and perhaps thermal conduction (fNar avan fc MedvedevI 
I2001t iChandran fc Maron|[20p: |Lazarian 200^). On an- 
other hand, in the outer region of a cluster, the thermo- 
dynamic state is likely to be dominated by gravitational 
processes and shock heating. Our models have included 
the essential physics in the outer regions. 

Figure [5] shows the ratio of electron and average ther- 
modynamic temperatures, r = T^/T, as a function of 
radius for the same cluster at z = 0. The shock ra- 
dius, i?sh ~ 4.2 Mpc, is about 1.5 times the virial ra- 
dius, i?vir ~ 2.8 Mpc. Within the virial radius, the 
temperature differences between electrons and ions are 
less than 1%, while beyond i?vir, T^/T decreases from 
~ 1 to 1/1800 at the shock radius for the model with 
P = 1/1800. For the /3 = 0.5 model, T^/f decreases from 
^ 1 to 0.5 at the shock radius. In general, Te/T decreases 
from ~ 1 to /3 at the shock radius for our models. Our 
models predict that Te/T ~ 0.8 at r « 0.9i?sh « 3.8 Mpc 
for a range of /3 between 1/1800 and 0.5. Beyond that 
radius, T^/T depends rather strongly on (3. The strong 
dependence at the shock radius can be used to distinguish 
shock heating models or constraint the shock heating ef- 
ficiency of electrons at the shock. 

The solid line in Figure |9] shows the projected 
spectroscopic-like temperature profiles, Tgi^non-cq, of the 
Msh = 1-53 X IO^^Mq cluster near the outer region at 
2 = for our non-equipartition model. We assume 
(3 — 1/1800 and Z = 0.3 Zq in this figure. For compar- 
ison, we also plot the projected spectroscopic-like tem- 
perature for a model where equipartition of electrons aiid 
ions throughout the cluster is assumed (Tg = Ti = T) 
as a dashed line. Both projected temperatures drop in 
the outer region as the radius increases, but the electron 
temperature in the non-equipartition model drops faster 
in most regions shown. Note that for the equipartition 
model, the projected temperature drop at the very last 
data point is a numerical artifact instead of a real fea- 
ture. This is due to the finite resolution of the shock 
handling in the hydro code which causes a slightly lower 
temperature compared to the idealized solution. Such an 
artifact does not affect our results significantly. 

The ratio of the projected temperature profiles of the 
two models, Tproj = Tsi,non-eq/7'si,cq, is plotted in Fig- 
ure [lOl Models for /3 = 0.5 and Z = 0.1 Zq are also 
plotted for comparison. Compared to the deprojected 
(or physical) temperature ratio (Figure ID), the deviation 
is larger for the projected temperature profiles, which is 
directly determined observationally. This is because for 
the projected temperature profile, electrons in the outer 
region also contribute to the inner region. For models 
with P = 1/1800 and Z = 0.3 Zq, there is a ~ 10% 
difference in the projected temperatures for the equipar- 
tition and the non-equipartition models at the virial ra- 
dius (i?vir ~ 2.8 Mpc) in contrast to less than a percent 
for the deprojected temperatures. The projected tem- 
perature difference increases to about 20% at a radius of 
~ 3.3 Mpc, which is about 1.2 of the virial radius. Even 
for /3 = 0.5, the difference at this radius can be as large 
as ~ 15%. We also found that the non-equipartition ef- 
fect on the projected temperature profiles is enhanced 
by metallicity. From Z = 0.1 to 0.3 Zq, the deviation 
at a radius of ^ 2.8(3.3) Mpc is enhanced by a factor of 
- 1.7(1.5) for the /3 = 1/1800 model. This is because the 



domination of the line emissions in the soft band spectra 
is enhanced by increasing metallicity, which is more im- 
portant for the non-equipartition model where electron 
temperature is lower. 

The effect of non-equipartition is larger for more mas- 
sive clusters. This is shown in Figure fTTl where Tpj-oj for 
clusters with different masses are plotted versus r/Rgh- 
This is because the deviation of the physical tempera- 
tures Te and T increases with cluster mass. The behavior 
qualitatively agrees with the analytic self-si milar solution 
in the Einstein-de Sitter universe found by iFox fc Loebl 
dHH). 

The surface brightness profiles for various energy bands 
of the Msh = 1-53 x 10^^ Mq cluster near the outer region 
at 2; = are shown in Figure fT2l The ratios S'non-cq/S'eq 
are plotted in Figure [131 We assume (3 = 1/1800 and 
Z — 0.3 Zq in these figures. For X-ray emission be- 
low ~ 2 keV, the relative difference between the non- 
equipartition and the equipartition models is smaller 
than that of the projected temperature profile. This 
is because the surface brightness depends on density 
squared but has a weaker dependence on temperature. 
The relative difference for X-ray emission above ^ 2 keV 
is similar to that of the projected temperature profile 
because of the energy cut off for lower temperature. 
The differences in surface brightness in all energy bands 
shown are < 10% for radii < 3 Mpc. Beyond ~ 3 Mpc 
where the electron temperature drops significantly be- 
low ~ 3 keV for the non-equipartition model, the dif- 
ferences in surface brightness profiles become important. 
The difference in X-ray surface brightness is most sig- 
nificant for the hard band (2.0-10.0 keV) compared to 
the soft (0.3-1.0 keV) and medium (1.0-2.0 keV) bands. 
For the soft band, the surface brightness for the non- 
equipartition model actually becomes larger than that of 
the equipartition model at large radii. This is because of 
the increase in the soft band emissivity for temperature 
below ~ 1 keV. We have also plotted the full X-ray band 
(0.3-10.0 keV) and the bolometric surface brightness pro- 
files in Figure [T2l The full X-ray band surface brightness 
profile for the non-equipartition model is always lower 
than that of the equipartition model, while the oppo- 
site is true for the bolometric surface brightness profile. 
This indicates that for the non-equipartition model, a 
large amount of emission occurs in the energy band be- 
low 0.3 keV compared to that of the equipartition model. 
In fact, the bolometric surface brightness near the shock 
radius for the non-equipartition model can reach ^ 35 
times that of the equipartition model. 

Current X-ray observations by Suzaku of cluster outer 
regions extend to only i?2oo ^ 2 Mpc where the non- 
equipartition effects on both the surface brightness and 
projected temperature profiles are < 1% in our mod- 
els. The sensitivities of these observations are limited 
by the Poisson variations in the background extragalac- 
tic sour ce density rather t han by the instrumental back- 
ground (jBautz et al.ll2009t ) . To push the sensitivity limit 
out to i?vir ~ 1 -4/^200 whcrc non-equipartition ef- 
fects on the projected temperature are ~ 10%, observa- 
tions combining high spatial resolution and high surface 
brightness sensitivity would be needed. On the other 
hand, increasing the coverage of solid angle can also 
help to reduce the Poisson variations. Here, we estimate 
how much improvement in sensitivity would be needed to 
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push the current limit of ^ i?200 out to i?vir ~ 1-4-R200- 
Assume most of the current ob servations only have sen- 
sitivities out to i?2oo (but see IGeoree et aLr[2009l ) and 
cluster surface brightness scales as shown in Figure [T2l 
From Figure [H 5'(i?2oo) ~ 6 5'(1.4i?2oo) for the 0.3- 
10.0 keV band, and hence to have a significant detection 
of X-ray emission at i?vir ~ 1.4i?200, a factor of ^ 6 
improvement in sensitivity would be needed. Probably 
a combination of the above two solutions (a factor of 6 
improvement in each case) is needed to push the sensi- 
tivity limit to ~ 1.4i?2oo- Pushing the limit to Rgh will 
be very chall enging with X-ray observations. Recently, 
iGeorge et al.l (j2009,,) reported that cluster emission has 
been detected out to ~ 1.5i?2oo in the cluster PKS0745- 
191. They have shown that the temperature at that ra- 
dius is ~ 30% lower than the temperature predicted by 
hydrodynamic simulations. If their results are confirmed, 
we suggest that this may be a signature of electron-ion 
non-equipartition. 

6.4. Results for the SZ Temperature Decrement 

Compared to X-ray observations, it is believed that 
cluster outer regions should be better studied by future 
SZ observations because the SZ effect depends only on 
the electron density to the first power, while X-ray emis- 
sion depends on the electron density squared. The X-ray 
signature drops much faster compared to the SZ signa- 
ture. Moreover, the SZ effect is independent of redshift, 
and a large number of high-redshift clusters should be 
observable. It has been suggested that future SZ ob- 
servations such as the Atacama Large Millimeter Array 
(ALMA) should be able to detect or rule out the pres- 
ence of accretion s hocks in clusters (jKocsis et al.l 120051 : 
iMolnar et all [20091 ). 

Figure [T31 shows the evolution of the temperature 
decrement magnitude due to the SZ effect, — ATsze, as 
a function of radius for a cluster with mass accreted to 
Msh = 1.53 X IO^^Mq at z = 0. The sohd line shows the 
non-equipartition model and the dotted line shows the 
equipartition model. Both the equipartition and non- 
equipartition models show very similar SZ temperature 
decrement profiles, with the SZ temperature decrement 
dropping faster in the non-equipartition models. For the 
cluster a.t z — 0, —ATsze drops from ^ 1 mK at about 
0.1 Mpc down to ^ 1 /xK at about 3.5 Mpc, and then 
drops very rapidly beyond that. For the same cluster at 
z = 2 when it had a mass of Mshiz = 2) = 6.72 x 
(corresponding to A/vir = Migg = 5.30 x at 
that redshift), the sharp drop in —ATsze occurs from 
a higher value of about 10 /uK near the shock radius. 
This suggests that the shock feature can be best studied 
through high-redshift clusters provided that the region 
of interest can be spatially resolved. To show the sig- 
nature of non-equipartition effect more clearly, we plot 
the difference between the SZ temperature decrements 
for the equipartition and the non-equipartition models, 
(5ATszE = ATszE,non-cq " ATszE,oq, in Figure [TSj It 
shows that the difference is larger at higher redshift for 
a given cluster. At z = 2, (5ATsze is of the order of 
1 fiK for the Msh(z = 2) = 6.72 x IO^^Mq cluster 
we considered. The ratios of the temperature decre- 
ments (or equivalently the Comptonization parameters), 
J/non-oq/j/oq_= ATszE,non-oq/ ATszE,oq , are also plotted 
in Figure 1161 Similarly to the temperature profiles, for 



z = at the virial radius of i?vir ~ 2.8 Mpc, Uaon-eq/ Veq 
is about 0.93, while it drops to about 0.8 at a radius of 
~ 3.5 Mpc. 

The detailed analysis of whether the effect can be 
distinguished observationally involves a discussion of 
the detail characteristic of the potential radio obser- 
vations, which will be given in an upcoming paper. 
Here, we only estimate roughly the possibility of whether 
the non-equipartition signature can be detected at the 
shock reg ion. W e closely f ollow the esti mation done 
by KocsislOn ([2005) and IMolnar et all §009). In 
their work, they studied whether two different models 
with and without a shock can be distinguished in fu- 
ture ALMA observations from the SZ effect. In prin- 
ciple, we can adopt their technique directly to distin- 
guish our two models with clusters in eq uipartition or 
not. In Figure 10(c) of IMolnar et al.l (|2009). they have 
argued that for two cluster models with JATsze of the 
order of a few /xK, the signal-to-noise ratio (S/N) for dis- 
tinguishing between their models can be as high as 70. 
For our model with the M^hiz = 0) = 1.53 x IO^^Mq 
cluster, ^ATszE is less than ~ 2 /xK even at z = 2. 
However, the most massive cluster today can reach a 
mass as high as A'hooiz = 0) w 2 x lO-^^M©. In Fig- 
ure [T71 we show the —ATsze profiles for the equiparti- 
tion and the non-equipartition models of a cluster with 
M2oo{z = 0) = 1.83 X lO^^M©, which corresponds to 
Mshiz = 0) = 3.06 X IO^^Mq. The profiles are simi- 
lar to those of the Msh{z = 0) = 1.53 x lO^^M© clus- 
ter, except the magnitudes are larger. The shock radius 
at z = 2 is about 1.5 Mpc, which corresponds to an 
angular size of 177" for the ACDM cosmology we as- 
sumed. The mass within the shock radius at z = 2 
is equal to Msh(z = 2) = 1.34 x IO^^Mq. Figure [111 
shows the deviation in ATsze between the equiparti- 
tion and non-equipartition models. At z = 2, the de- 
viation is i5ATsz E ~ 4 — 5 AtK near the shock radius. 
Closel y following iKocsis et al.l (|2005D and IMolnar et al.l 
([2009'), we estimate the S/N for distinguishing between 
the equipartition and the non-equipartition models of the 
very m assive cluster at z = 2 by ALMA. IKocsis et al.l 
(|2005f) estimated that a 20 hr on-source integration 
time for ALMA would be enough to achieve a sensitivity 
of < 10/iK, but the latest ALMA Sensitivity CalculatoiQ 
shows that about 260 hr would be needed for the required 
angular resolution (FWHM) of 2" at 100 GHz with 64 an- 
tennas. Though such a very long observation is possible, 
we lower the required angular resolution (FWHM) to 4" 
to increase the sensitivity, and the required on-source in- 
tegration time is reduced to ~ 17 hr. Hence we assume a 
~ 17 hr on-source integration time and estimate the S/N 

1/2 

ratio to be « N^l^ (S/N)i, where iVpix is the number of 
independent pixels in the region of interest, and (S/N)i is 
the S/N for a single pixel in ALMA. From Figure [TH the 
width of the annular region of interest from ^ 1 to ^ 1.5 
Mpc is about 0.5 Mpc ~ 59" at z = 2, where 1 Mpc 
corresponds to 118" at this redshift for the standard 
ACDM universe assumed. Assuming a Gaussian beam 
with an angular diameter (FWHM) of 4" for ALMA, the 
number of beams within the region of interest are about 
iVpix ~ 3000. In a real cluster, the accretion would likely 



i:/ /www. eso.org/sci/facilities/alma/observing/tools/etc/index. htm] 



14 



Wong & Sarazin 



not to be spherical, and material would accrete through 
filaments. However, in regions other than the filaments, 
accretion shocks are roughly spherical and those are the 
region s which our model may be applied. iMolnar et al.l 
([200!^ have estimated the area coverage factor of such 
spherical accretion shock region to be ~ 50%. With 
this correction, iVpix,5o% ~ 1500. The noise for a sin- 
gle pixel is estimated to be 10/zK (jKocsis et al.l 120051 : 
IMolnar et~a l. 2009). We take the signal for a single pixel 
to be Si sa 5ATsze ~ 4^K (Figure [T5)) . Thus, we obtain 
an S/N of ^ 16. Even if the area coverage factor goes 
down to only 10%, the S/N c an still be as high as 7. As 
noted in IMolnar et al.l (|2009l ) , the sensitivity of interfer- 
ometers is reduced for large-scale smooth surface density 
distributions, which is perhaps the case for the SZ effect 
in the cluster outer regions. The S/N in this analysis may 
be overestimated. However, they also argue that such a 
reduction in sensitivity can be recovered by using non- 
linear de-convolution of data from mo saic observations 
(jHelfer et al.ll2002l : IMolnar et al.ll2009f ). The above esti- 
mation assumes that the cluster model parameters (e.g., 
cluster mass, shock radius) are known in advance. In 
real observations trying to distinguish between equipar- 
tition and non-equipartition models, if those parameters 
cannot be obtained by other means, we may need to fit 
the parameters from the data. This will in general re- 
duce the S /N estimated above. We defer a more detailed 
study to a future paper, but based on the rather high 
S/N estimated from the arguments above, we conclude 
here that future SZ observations, such as those done by 
ALMA, may be able to distinguish between equipartition 
and non-equipartition models near the shock region. 

6.5. Results for the Integrated SZ Biases, Ynon-cq/Ycq, 
and its Evolution 

Figure ITOl shows the integrated SZ biases, Fnon-oq/^q 
defined in Section 16.2! as a function of Mgh for both our 
simulated realistic NFW model in the ACDM universe 
and the numerical simulated self-similar model in the 
Einstein~de Sitter universe at different redshifts. We as- 
sume /gas = 0.17 for models in the Einstein-de Sitter 
universe in Figure [T21 The masses are evaluated at the 
labeled redshift so that the evolutionary history for a 
particular cluster cannot be seen on the graph. Both 
models show that ynon-cq/X^q increases as the cluster 
mass increases, which is expected as the effect of non- 
equipartition increases with mass. For the realistic NFW 
model, ynon-eq/i^q dccrcascs from 1 for Mgh = IO^^Mq 
down to - 0.9 for M^h = IO^^Mq at z = 0. The up- 
per mass limit to which a cluster can grow is limited by 
the background cosmology as well as the initial density 
fluctuation amplitudes. The most massive nearby clus- 
ter observed has a mass of about M200 ~ 2 x 10^^ Mq, 
which corresponds to M^h » 3.3 x lO^^M© in our NFW 
model. Ynon-cq/yeq for the most massive cluster in our 
universe is hence about 0.97 if it is nearby (z « 0). This 
bias would be larger (ynon-oq/^oq smaller in magnitude) 
at higher redshift for a given mass. 

Recent observations suggest that a significant fraction 
(20% -40%) of the thermal energy is m i ssing from clus- 
ters (lEttoril [20 03: LaRoque '"eral1 120061: lyikhlinin et al.l 
[20061 lAfshordi et al. 2q07|; lEvrard et al.l 120081 but lee 
also [Giodini et al.. .200^x 7 Obviously, if electrons and 



ions are in non-equipartition, the thermal energy mea- 
sured by X-ray or SZ observations should be reduced. 
Our simulations suggest that for cluster with Mgh ~ 
1.5 X lOi'''M0(M2oo IO^^Mq), the non-equipartition 
effect can account for only about 2%-3% of the missing 
thermal energy. For the most massive clusters, up to 3%- 
4% of the thermal energy beyond the equipartition value 
may be stored in the ions near the shock radius, if elec- 
trons and ions are in non-equipartition. For Msh smaller 
than about 5 x at z = 0-2, the non-equipartition 

effect is less than 1%. 

For the self-similar model, the ynon-oq/ioq curves at 
the four different redshifts actually lie almost along the 
same line which cannot be easily distinguished from the 
graph (Figure [T9l) . This shows that the integrated SZ 
bias for the self-similar model does not evolve with red- 
shift. In contrast, the integrated SZ bias for the real- 
istic NFW model in the ACDM universe evolves with 
redshift. In the ACDM universe, the expansion of the 
universe starts to accelerate around the redshift where 
flM ^ ^A', this breaks the self-similar solution for cosmo- 
logical accretion, and hence we should expect ynon-oq/^^q 
would also deviate from self-similarity in general. For 
our realistic NFW model in the ACDM universe, the in- 
tegrated SZ bias decreases as z decreases. This is proba- 
bly due to the decreasing rate of accretion onto clusters 
in the ACDM universe during the cosmological acceler- 
ation, which results in a relatively longer timescale for 
the electron-ion equilibration inside a cluster compared 
to a cluster with the same mass in the Einstein-de Sit- 
ter universe. At z = 0, Fnon-eq/^eq of the realistic NFW 
model is smaller than that of the self-similar model with 
the same /gas for Msh > IO^^Mq, but the effect is similar 
for Msh < lO^^M© for both models. At higher redshifts, 
the effect of non-equipartition is larger for the realistic 
NFW model for the entire mass range. 

Though the magnitude of Fnon-cq/icq is small for the 
range of cluster masses, even a percentage level deviation 
in the most massive clusters is important for precision 
cosmology studies. Cosmological studies using the mass 
function evolution depend sensitively on the high mass 
end of clusters. For cluster mass count surveys using Y 
as a mass proxy, if the bias in Y is not properly taken into 
account, the resulted mass function would be biased low 
at the high mass end (i.e., less massive clusters would 
be observed if clusters are in non-equipartition). Even 
though the mass-F relation can be self-calibrated, the 
evolution of Ynon-cq/Yeq may introduce a bias if the self- 
calibration is not properly done at each redshift. For 
example, if the mass and Y relation is self-calibrated 
correctly for low-redshift clusters but this calibration is 
extrapolated to high-redshift clusters, bias would be in- 
troduced if the non-equipartition effect is not properly 
taken into account. A detailed study of the effect of 
non-equipartition on SZ surveys and the implication to 
cosmological studies will be presented in an upcoming 
paper. 

7. DISCUSSIONS AND CONCLUSIONS 

Using one-dimensional hydrodynamic simulations, we 
have calculated a sample of realistic NFW clusters in 
a range of masses in the ACDM cosmology. The clus- 
ter properties we simulated are consistent with the one- 
dimensional A^-body simulations bv .Rvu &: Kane. (,1997.) . 
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and they have shown that their calculations repro- 
duce the density and temperature profiles of the three- 
dimensional simulated relaxed clusters in the outer re- 
gions. Our one-dimensional hydrodynamic simulations 
help us to isolate the important physical processes under 
controlled conditions. We have studied in detail the ef- 
fect of non-equipartition in the outer regions of relaxed 
clusters in the ACDM cosmology. 

Using /gas = 0.17 (which is the upper limit for a 
cluster), we give a conservative lower limit of the non- 
equipartition effect on clusters. We have shown that for 
a cluster with a mass of Afsh ^ 1.5x lO^^M©, within i?vir, 
electron and ion temperatures only differ by less than a 
percent. Our results show that the effect is smaller than 
those calculated from recent three-dimensional simula- 
tions, which shows that Tg can be biased low by 5% at 
i?2oo ~ 0.7i?vir (model CL104 in Rudd & Nagai 2009). A 
detailed analysis is needed to address the difference, but a 
possible explanation may be that in a three-dimensional 
cluster the accretion shock can be formed further in. Our 
results show that Tg/T can reach « 0.8 for a range of 
non-adiabatic electron heating efficiency /? ~ 1/1800 to 
0.5 at ~ 0.9i?sh (or ~ 1.4i?vir)- Beyond that radius, 
Tg/T depends rather strongly on /3, and such a strong 
dependence at the shock radius can be used to distin- 
guish shock heating models or constraint the shock heat- 
ing efficiency of electron. We also show that the effect 
of non-equipartition is larger for more massive clusters, 
which is consistent to analy tic self-similar mod els in the 
Einstein~de Sitter universe (iFox fc LoeblfTOOl. 

Using the algorithm developed bv IVikhlininI (|2OO60 
which takes into account the soft emission at low tem- 
perature down to ~ 0.5 keV, arbitrary metallicity, 
and instrumentation response, we calculated the X-ray 
spectroscopic-like temperature profiles which are the one 
to be directly determined observationally. The effect of 
non-equipartition on the projected temperature profiles 
is larger than that on the deprojected (or physical) tem- 
perature profiles. Non-equipartition effects can introduce 
a ~ 10% bias in the projected temperature at i?vir for a 
wide range of /?. This is because for the projected tem- 
perature profile, electrons in the outer region also con- 
tribute to the inner region. We also found that the effect 
of non-equipartition on the projected temperature pro- 
files can be enhanced by increasing metallicity. This is 
because the domination of the line emissions in the soft 
band spectra is enhanced by increased metallicity, which 
is more important for the non-equipartition model where 
the electron temperature is lower. 

The effect of non-equipartition on X-ray surface bright- 
ness profile in the 0.3-2 keV band is smaller than that 
on the projected temperature profile. This is because the 
surface brightness depends on density squared but with 
a weaker dependence on temperature. This means that 
in the outer regions, clusters in non-equipartition have 
similar X-ray surface brightness profile for i? < 2 keV, 
but with bigger difference in temperature compared to 
those equipartition counterparts. For E > 2 keV, non- 
equipartition effects on X-ray surface brightness profiles 
are similar to those on the projected temperature pro- 
files. For a cluster with M^h ~ 1.5 x IO^^Mq, the effect 
of non-equipartition on surface brightness profiles in all 
energy bands is < 10% for radii < 3 Mpc; beyond that, 



the effect can be important. We found that for the non- 
equipartition model, the surface brightness profile in the 
low-energy band < 1 keV can be higher than that of the 
equipartition model in the cluster outer regions. Non- 
equilibrium ionization, which was not considered in our 
calculations of the emissivities, may further enhance the 
line emissions in the soft bands {E <1 keV). 

Current X-ray observations extend to only ~ i?200 ~ 
2 Mpc, although some results from recen t Suzaku ob- 
scrvations begin to go a bit beyond that ([George et al.l 
2009). Within those regions with < i?2oo, electrons and 
ions should be almost in equipartition and the signatures 
in the X-ray temperature and surface brightness should 
be rather weak. But future X-ray observations may ex- 
tend to ~ i?vir ~ l-4i?200 or cvcu close to the shock ra- 
dius. We have shown that non-equipartition of electrons 
and ion s should be de tectable in those studies. The re- 
sults by iGeorge et al.l ([2009) support our conclusion. 

The effects of non-equipartition on the SZ effect were 
studied. At z = 0, the effect on the Comptonization 
parameters is similar to that of the projected tempera- 
ture profiles. For a cluster with Al^h ^ 1.5 x IO^^Mq, 
J/non-cq/j/cq ~ 0.93(0.8) at l(1.3)i?vir- For a given cluster, 
the difference between the SZ temperature decrements 
for the equipartition and the non-equipartition models 
is larger at a higher redshift. For the most massive 
clusters at z w 2, the differences can be SATsze ~ 4- 
5 /iK near the shock radius. A detailed analysis of 
whether the equipartition and non-equipartition models 
near the shock region can be distinguished by, for exam- 
ple, ALMA, will be presented in a future paper. 

The effects on the integrated SZ Comptonization pa- 
rameter, which measures the thermal energy content of 
the electrons, were studied. We have shown that the 
integrated SZ bias, ynon-oq/^cq, increases as the cluster 
mass increases, which is expected as the effect of non- 
equipartition increases with mass. In general, the non- 
equipartition effect is larger for the realistic NFW model 
in the ACDM universe than that for the self-similar 
model in the Einstein-de Sitter universe, assuming that 
they have the same /gas- Our simulations suggest that 
for relaxed clusters with Mgh ~ 1.5 x IO^^Mq, the non- 
equipartition effect can account for only about 2%-3% of 
the missing thermal energy globally. For the most mas- 
sive clusters, up to 4%-5% of the thermal energy beyond 
the equipartition value may be stored in the thermal en- 
ergy of ions near the shock radius, but for clusters with 
^^sh ^ 5 X 1O^*M0, the non-equipartition effect is less 
than 1%. Thus, we argue that, at least for relaxed clus- 
ters, the non-equipartition effect alone can only account 
for some of the missing thermal energy problem, if any, 
for high-mass clusters but not for clusters with smaller 
masses. On the other hand, this suggests that hot gas 
may be missing due to other astrophysical processes not 
yet known, and the /gas in a real cluster should be lower 
than that we used in our numerical simulations. We have 
estimated that reducing /gas by 20% will enhance the lo- 
cal non-equipartition effect near the outer region by a 
few percent, but the integrated SZ bias, Fnon-cq/^cq is 
not affected by more than a percent. 

We emphasis here that even though non-equipartition 
effects may not affect the global energy budget signifi- 
cantly, the effect is still important locally in the outer 
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regions (~ i?vir) of a cluster. Future X-ray and SZ ob- 
servations may extend out to i?sh, and the effect of non- 
equipartition should be considered when studying cluster 
properties in those regions. 

We found that for our realistic NFW model in the 
ACDM universe, i^on-eq/i^oq evolves with redshift, which 
is in contrast to the self-similar model in the Einstein- 
de Sitter universe. For our realistic NFW model in the 
ACDM universe, Ynon-cq/Yeq decreases as z decreases. 
This is probably due to the decreasing rate of accretion 
onto clusters in the ACDM universe during the period 
of cosmological acceleration, which results in a relatively 
longer timescale for the electron-ion equilibration inside 
a cluster compared to a cluster with the same mass in 
the Einstein-de Sitter universe. Though the magnitude 
of Ynan-eq/Yeq is Small for the range of cluster masses, 
even a percentage level deviation in the most massive 
clusters can be important for precision cosmology stud- 
ies. Such a variation of Y with z would introduce an 
apparent evolution in /gas, which would bias the cos- 



mological studie s using the /gas techn iques (| Allen et al.l 
I2008f) . Recently, iRudd fc Nagal ()2009D have shown that 
the non-equipartition effect on Y can be enhanced by ma- 
jor mergers up to 30%, although for low Mach number 
mergers, the shock heating efficiency for electrons may 
be higher which can weaken the non-equipartition effect 
(|Ghavamian et al.l 120071: iMarkevitch fc Vikhlinin| [2007'). 
The temporary boost due to mergers may have a signifi- 
cant effect on the estimation of cosmological parameters 
using clusters. We defer a detailed study of the effect on 
cosmology studies in a future paper. 
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APPENDIX 

EXPRESSION FOR ELECTRON HEATING WITHIN THE SHOCK 
Adiabatic changes in electron temperature are given by 

(T'e2)ad = ( ) ^el > 



(Al) 



where Pe2/Pei — Pg2/Pgi for fully ionized plasma. In addition to the adiabatic heating, electrons can also be heated by 
non-adiabatic processes. From the definition of non-adiabatic electron heating efficiency (equation (|26p ). the change 
of electron temperature due to non-adiabatic heating is 



non-ad — 

Thus, the final electron temperature can be expressed as 

Te2 = Tel + (ATe)non-ad ■ 

\PgiJ 

Similarly for the average thermodynamic temperature, we have 



(T2)ad = 



Pg2 
Pgl 
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Combining equations (|A2p . (|A3p . and (|A5p . we get 

T,.2=(^) Tei+/3max 

31/ 



0,T2 



Pg2 
Pgl 
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Ti 



(A5) 
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where a minimum of zero in the second term is set to ensure that numerical fiuctuations do not introduce false decreases 
in entropy in the non-adiabatic heating. 
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Fig. 1. Gas velocity (upper panel), density (middle panel), and pressure (lower panel) profiles of our simulated cluster for the self- 
similar coUisionless dark matter dominated accretion model in the Einstein-de Sitter universe at four different rcdshifts. The model has 
a dark matter mass accreted within Rirs of AIns = IO^^Mq at z = 0. The insets in the upper and middle panels show the large radius 
behavior of the gas, with the velocity profile shown for z = only. The circles on the inset of the velocity profile give the Hubble flow 
velocity at z = 0. We set Qf, = 0.05 in this model. 
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Fig. 2. — Scaled gas velocity (upper panel), density (middle panel), and pressure (lower panel) profiles of our simulated cluster. The 
hydrodynamic model is the same as in Figure^ All four lines lie almost on top of one another in the figures. Circles show the corresponding 
analytic solutions. 
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Fig. 3. — Similar to Figure[2]but for the self-similar coUisional gas dominated accretion model in the Einstein-de Sitter universe. The 
model has a gas mass accreted within Rirs of Afi78 = 1.18 X 10^^ A/0 at 2 = 0. 



Non-Equipartition Signatures in Cluster Outskirts 



21 




Fig. 4. — Ratio of the electron and average thermal dynamic temperature, r, as a function of scaled radius r/i?sh for clusters of the 
self-siniilar coUisionless dark matter dominated accretion model in the Einstein-de Sitter universe at 2; = 0. The upper and lower solid 
lines are the analytic solutions for cluster masses with Miyg = 5.26 X and Mi7g = 1.05 X respectively. The open circles 

on the corresponding lines are our simulated results. We assume Cli, = 0.05 here. 
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Fig. 5. — Gas velocity (upper panel), density (middle panel), and pressure (lower panel) profiles of our simulated cluster for the realistic 
NFW-dark energy model in the standard ACDM universe at four different rcdshifts. The model has total mass accreted within iigij of 
A^sh = 1-53 X 10^^ Mq at z = 0. The insets show the large radius behavior of the gas, with the velocity profile shown for z = only. The 
circles on the inset of the velocity profile give the Hubble flow velocity at z = 0. 
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Fig. 6. — Average thermodynamic temperature profiles of our simulated realistic NFW-dark energy model cluster. The hydrodynamic 
model is the same as in Figure [3] 
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Fig. 7. — Average thermodynamic temperature profiles of our simulated self-similar dark matter dominated cluster. The hydrodynamic 
model is the same as in Figure [Tl 




Fig. 8. — Ratio of electron and average thermodynamic temperatures, t = Te/T, as a function of radius at z = 0. The hydrodynamic 
model is the same as in Figure [5] The models with shock heating efficiency /3 = 1/1800 and 0.5 are shown in solid and dashed lines, 
respectively. 
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Fig. 9. — Projected X-ray spectroscopic-like temperature profiles for cluster models at z = 0. The hydrodynamic model is the same as 
m Figure [5] The non-equipartition model with /3 = 1/1800 is shown in solid line, while the equipartition model is shown in dashed line. 
We assume Z = 0.3 Zq here. 
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Fig. 10. — Ratio of the projected X-ray spectroscopic-like temperature profiles of the non-equipartition (with j3 = 1/1800 and 0.5) and 
the equipartition models at z = 0. Models with Z = 0.1 and 0.3 Zq are presented. The hydrodynamic model is the same as in Figure[5] 
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Fig. 11. — Ratio of the projected X-ray spectroscopic-like temperature of the non-equipartition with /3 = 1/1800 and the equipartition 
models as a function of r/Hsh at z = for different cluster masses. The hydrodynamic models are similar to those in Figure [S] but with 
Msh = 7.65, 15.3, and 30.6 x IO^^Mq at z = 0. We assume Z = 0.3 Zq here. 
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Fig. 12. — Surface brightness profiles for various energy bands for tlie non-equipartition (solid lines) and equipartition (dashed lines) 
cluster models at ^ = 0. Models are the same as in FigurelQ] 
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Fig. 13. — Ratios Snon-eq/Scq as a function of radius. Models are the same as in Figurc [T2l The ratio of the bolometric surface brightness 
near the shock radius reaches ~ 35 (outside the scale of the figure). 
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Fig. 14. — SZ temperature decrement profiles, — ATgzEi of our simulated cluster at four different redshifts. The non-equipartition models 
(/3 = 1/1800) are shown in solid lines, while the equipartition models are shown in dashed lines. The hydrodynamic model is the same as 
in Figure [5] 
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Fig. 15. — Difference (5ATsze between the SZ temperature dec rements of the equipartition and the non-equipartition models (/3 = 1/1800) 
at four different redshifts. Models are the same as in Figure [141 




Fig. 16. — Ratio ?/non-cq/?/cq between the Comptonization para meters of the equipartition and the non-equipartition models (/3 = 1/1800) 
at four different redshifts. Models are the same as in Figure lT4l 




r (Mpc) 



Fig. 17.— Same as Figure [Til but with M^hC^ = 0) = 3.06 x W^^Mq for the cluster model. 
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Fig. 18.— Same as Figure [TSl but with M^hC^ = 0) = 3.06 X IO^^Mq for the cluster model. 




Fig. 19. — Integrated SZ biases, Viion-cq/^^cq, as a function of Msh for both our simulated realistic NFW model in the ACDM universe 
(thin lines) and the numerical simulated self-similar model in the Einstein-de Sitter universe (thick lines) at four different redshifts. We 
assume /gas = 0.17 for models in the Einstein-de Sitter universe. The four lines for the self-similar model lie almost along the same line 
which cannot be easily distinguished on the graph. 
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TABLE 1 

Masses and Radii'' of Some Representative NFW 
Cluster Models in the Standard ACDM Cosmology 

AT 2 = 







vir 


-Rvir 


M2()0 


-R2()() 




-R.soo 


7.65 
15.3 
30.6 


3.31 
4.22 
5.41 


6.05 
11.9 
23.1 


2.19 

2.75 
3.43 


4.90 
9.50 
18.3 


1.60 
1.99 
2.48 


3.59 
6.83 
12.9 


1.06 
1.32 
1.63 



Masses are in the unit of 10 Mq and radii are in 
the unit of Mpc. 

Avir = 95.3 at z = for the standard ACDM cos- 
mology. 



